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Global Existence and Large Time Asymptotic Behavior of Strong 
Solutions to the Cauchy Problem of 2D Density-Dependent 
Navier-Stokes Equations with Vacuum* 

Boqiang Lii^ Xiaoding Shi^ Xin Zhong^ 


Abstract 

We are concerned with the Cauchy problem of the two-dimensional (2D) nonhomogeneous 
incompressible Navier-Stokes equations with vacuum as far-held density. It is proved that 
if the initial density decays not too slow at inhnity, the 2D Cauchy problem of the density- 
dependent Navier-Stokes equations on the whole space admits a unique global strong 
solution. Note that the initial data can be arbitrarily large and the initial density can contain 
vacuum states and even have compact support. Furthermore, we also obtain the large time 
decay rates of the spatial gradients of the velocity and the pressure which are the same as 
those of the homogeneous case. 
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1 Introduction 

The motion of a two-dimensional nonhomogeneous incompressible fluid is governed by the fol¬ 
lowing Navier-Stokes equations: 

dtp -I- div(/?n) = 0, 

< dt{pu)+ diw{pu®u) — pI^u + VP = (1.1) 

divu = 0. 

\ 

Here, t > 0 is time, x G is the spatial coordinates, and the unknown functions p = p{x, t), u = 
{u^ and P = P{x,t) denote the density, velocity, and pressure of the fluid, respectively; 
p > 0 stands for the viscosity constant. 

We consider the Cauchy problem of (|l.ll) with {p, u) vanishing at infinity (in some weak sense) 
and the initial conditions: 


p{x,Q) = Pq{x), pn(x,0) = PqUq{x), x G M^, 
for given initial data po and uq. 


( 1 . 2 ) 
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The study of the system dni) has a long history. Since the pioneering work of Leray [Ml, 
there are huge literatures on the studies of the large time existence and behavior of solutions of 
(HH) in the case that the density is constant (refer to [5l[T3l[T5l[T9l[25] and references therein). 
The mathematical study of nonhomogeneous incompressible flow goes back to the seventies of the 
last century. When the initial density is strictly away from vacuum, Kazhikov [13] established 
the global existence of weak solutions (see also [H). Later, Antontsev-Kazhikov-Monakhov [2| 
gave the first result on local existence and uniqueness of strong solutions. Moreover, they also 
proved that the unique local strong solution is a global one in two dimensions. Recently, the 
global existence result of dni) with small initial data belonging to certain scale invariant spaces 
was obtained in [7j. 

On the other hand, when the initial density allows vacuum, Simon [23| proved the global ex¬ 
istence of weak solutions, which was extended by Lions [19] to the case of density-dependent vis¬ 
cosity. For the initial density allowing vacuum, Choe-Kim [3] proposed a compatibility condition 
and obtained the local existence of strong solutions of (dH) with large data in two-dimensional 
bounded domains. Recently, Huang-Wang m showed that the local strong solution obtained 
in [3] is indeed a global one. Very recently, motivated by [laisi], Lii-Xu-Zhong [20] establihsed 
the local (with generally large data) existence of strong solutions to 2D Cauchy problem of the 
nonhomogeneous incompressible magnetohydrodynamic equations with vacuum as far-held den¬ 
sity, which in particular yields the local existence of strong solutions to the 2D Cauchy problem 
of (11.11) . However, the global existence of strong solution to the 2D Cauchy problem of (11.11) with 
vacuum and general initial data is still open. In fact, this is the main aim in this paper. 

Before stating the main results, we hrst explain the notations and conventions used throughout 
this paper. For i? > 0, set 

Br = {x G |x| < R} , / -dx = / -dx. 

J Jm? 

Moreover, for 1 < r < oo and /c > 1, the standard Lebesgue and Sobolev spaces are dehned as 
follows: 






Rk ^ 


Now we dehne precisely what we mean by strong solutions. 

Definition 1.1 If all derivatives involved in dni) for (p, u, P) are regular distributions, and 
equations dni) hold almost everywhere in x (0,T), then {p,u,P) is called a strong solution to 

CU). 


Without loss of generality, we assume that the initial density pQ satishes 



= 1 , 


which implies that there exists a positive constant Nq such that 


(1.3) 


Our main result can be 

Theorem 1.1 In addition 
given numbers a > 1 and q 

Po > 0, pox°- 



stated as follows: 

to ()1.3p and ()1.4I) . assume that the initial 


> 2 , 

GL^nH^n divuo = 0, Vuo G L^ 


(1.4) 


data {po,uo) satisfy for any 
y/f^uo G L^, (1.5) 


2 


where 


X = {e + \x\'^)^^‘^\og^{e + \x\^). ( 1 . 6 ) 

Then the problem (inii-dol) has a unique global strong solution (p, u, P) satisfying that for any 
0 < T < oo, 

'0<p£ C([0, T];L^nH^n 
px^ e L°°(0, T; 

yfpu, Vu, x~^u, \/i^ut, y/iVP, \/tV^u G L°°(0, T; L^), 

< Vu £ L‘^{0,T;H^) n ( 1 . 7 ) 

VP G L2(0, T; L2) n L(''+i)/'?(0, T; L1), 

VtVu G L2(0,T;1V^’''), 

VtVut, y/tx~^ut G L^(M^ x (0,r)), 

and 

inf / p(x,t)dx>-, (1-8) 

/or some positive constant Ni depending only on \\y/f^uo\\i 2 ,No, and T. Moreover, {p,u,P) has 
the following decay rates, that is, for t >1, 

\\Vu{;t)\\l, + \\vM-Mp + \NP{-Mp < ct-\ ( 1 . 9 ) 

where C depends only on p, ||/9o||LinL°°) II\/)^^o||l2; o-n^d ||Vno||L2- 

Remark 1.1 Our Theorem \l.l\ holds for arbitrarily large initial data which is in sharp contrast 
to Li-Xin m where the smallness conditions on the initial energy is needed in order to obtain 
the global existence of strong solutions to the 2D compressible Navier-Stokes equations. 

Remark 1.2 It should be noted here that although the equations m2 degenerate near vacuum, 
our large time decay rates (USD are the same as those of the homogeneous case In- 

Remark 1.3 Compared with fll]/ . there is no need to impose the additional compatibility condi¬ 
tions on the initial data for the global existence of the strong solution. 

We now make some comments on the key ingredients of the analysis in this paper. Note that 
for initial data satisfying m, Lii-Xu-Zhong [20] established the local existence and uniqueness 
of strong solutions to the Cauchy problem (ll.ll) - (ll.2ll (see Lemma [XT]) . Thus, to extend the local 
strong solution to be a global one, one needs to obtain global a priori estimates on strong solutions 
to (ll.ip - (ll.2l) in suitable higher norms. It should be pointed out that the crucial techniques of 
proofs in m cannot be adapted directly to the situation treated here, since their arguments 
depend crucially on the boundedness of the domains. Moreover, it seems difficult to bound the 
L'^(M^)-norm of u just in terms of || v^^IIl 2 (r 2 ) and ||Vri|| 2 , 2 (R 2 ). To overcome these difficulties 
mentioned above, some new ideas are needed. First, motivated by [HlUllIH]; multiplying (|l.ll) 9 
by the material derivatives of the velocity, it = ut + u ■ Vu, instead of the usual ut (see m) , we 
find that the key point to obtain the estimate on the L°°(0, T; L^(M^))-norm of the gradient of 
the velocity is to bound the term 




PdiU^ djU^dx. 


Motivated by |6|, we hnd that I 2 in fact can be bounded by ||VP||^21| Vri ||^2 (see (]3.8l) i since 
diMdjU^ G due to the facts that divu = 0 and that V-*- • Vu = 0 (see Lemma [23] also). Next, 
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to obtain the estimates on the gradient of the density, motivated by [H HlTTlfT^ . we apply the 
operator dt + u-V to (jl.ip n and multiply the resultant equality by ii^ to get the time-independent 
estimates on both the L°°(0, T ; L^(M^))-norm of and the L^(R^ x (0, T))-norm of 

(see ()3.13p l which combined with some careful analysis on the spatial weighted estimate of the 
density (see (13.211) 1 thus yield the bound on the L^(0,T; L°°)-norm of the gradient of the velocity 
(see (13. 32 p i which in particular implies the desired bound on the L°°{0, T ; L'^)-norm of the gradient 
of the density. Finally, with these a priori estimates on the gradients of the density and the velocity 
at hand, one can estimate the higher order derivatives by using the same arguments as in [lllfT^ 
to obtain the desired results. 

The rest of this paper is organized as follows. In Section [2j we collect some elementary facts 
and inequalities that will be used later. Section [3] is devoted to the a priori estimates. Finally, 
we will give the proof of Theorem 11.11 in Section 01 

2 Preliminaries 

In this section we shall enumerate some auxiliary lemmas used in this paper. 

We start with the local existence of strong solutions whose proof can be found in [201 Theorem 

1 . 2 ]. 

Lemma 2.1 Assume that {pq,uq) satisfies (jl.Sp . Then there exist a small time T > 0 and a 
unique strong solution {p,u,P) to the problem (|1.1I) - (I1.2I) in X (0, T) satisfying (11.71) and 

im- 

Next, the following Gagliardo-Nirenberg inequality (see [22]) will be used later. 

Lemma 2.2 (Gagliardo-Nirenberg) For q G [2,oo),r G (2,oo), and s G (l,oo), there exists 
some generic constant C > 0 which may depend on q, r, and s such that for f G and 

g G L"(M2) n Dh''(R2), we have 

ll/ll?,,,*.) 

The following weighted L”* bounds for elements of the Hilbert space Z)^’^(R^) = {u G 
G L^(R^)} can be found in [T9l Theorem B.lj. 

Lemma 2.3 For m G [2,oo) and 0 G (1 -|- m/2,oo), there exists a positive constant C such that 
for all V G ii^’^(R^), 

) l/m 

< CIHlFBi) + C'l|Vu||i2(R2). (2.1) 

The combination of Lemma [2.3l and the Poincare inequality yields the following useful results 
on weighted bounds, whose proof can be found in m Lemma 2.4]. 

Lemma 2.4 Let x be as in ([ISD. Assume that p G L^(R^) n L°°(R^) is a non-negative function 
such that 



for positive constants Mi, M 2 , and Ni > 1. Then for e > 0 and g > 0, there is a positive constant 
C depending only on e,77, Mi,M2, and Ni, such that every v G L>^’^(R^) satisfies 

^llL(2+e)/^)(R2) ^ -I- C11 Vu11i;^2(r2), (2.2) 

with fj = min{l, g}. 
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Finally, let and BMO(M^) stand for the usual Hardy and BMO spaces (see [JU Chapter 

IV]). Then the following well-known facts play a key role in the proof of Lemma 13.21 in the next 
section. 


Lemma 2.5 (a) There is a positive constant C such that 




for all E G and B G satisfying 

divF; = 0, V-^-H = 0 in 
(b) There is a positive constant C such that 

lbllBMO(IR2) < C||Vu||2,2(]r2), 

for all V G 

Proof, (a) For the detailed proof, see [U Theorem 11.Ij. 

(b) It follows from the Poincare inequality that for any ball H C 



which directly gives (12.3p . 



dx < C 



(2.3) 


□ 


3 A Priori Estimates 

In this section, we will establish some necessary a priori bounds for strong solutions (/?, u, P) to 
the Cauchy problem (|l.lD - (ll.2p to extend the local strong solution. Thus, let T > 0 be a fixed 
time and {p,u,P) be the strong solution to (ll.ip - (ll.2|) on x (0, T] with initial data {po,uo) 
satisfying (ll.3D - (ll.5p . 

In what follows, we will use the convention that C denotes a generic positive constant de¬ 
pending on /r, a, and the initial data, and use C{a) to emphasize that C depends on a. 


3.1 Lower Order Estimates 


First, since divu = 0, we state the following well-known estimate on the L°°(0, T; L^)-norm of 
the density. 

Lemma 3.1 ( m) There exists a positive constant C depending only on ||/Oo||Lini/°° such that 

sup IjpIlLinL- < C. (3.1) 

te[o,T] 

Next, the following lemma concerns the key time-independent estimates on the L°°{0,T; L'^)- 
norm of the gradient of the velocity. 

Lemma 3.2 There exists a positive constant C depending only on p, ||/?o|| 1^00 ^ ||v^uo||l 2 , and 
||Vuo||l 2 such that 


sup ||Vu ||^2 + / p\u\‘^dxdt < C, 

te[o,r] Jo J 

where ii = ut + u - Vu is the material derivatives of the velocity. Furthermore, one has 

rT 


sup t||Vri |||2 + [ t f p\u\'^dxdt < C. 
te[o,r] Jo J 


(3.2) 


(3.3) 
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(3.4) 


Proof. First, applying standard energy estimate to (inD gives 

sup \\y/pu \\\2 + f \\Vu\\l 2 dt < C. 
t£[0,T] Jo 

Next, multiplying m 2 by u and then integrating the resulting equality over lead to 

J p\u\‘^dx = j{pAu ■ ii — VP • u)dx = Ii + 12 - (3.5) 

It follows from integration by parts and Garliardo-Nirenberg inequality that 

Ii = J pAu ■ {dtu + u ■ Vu)dx 

=-^^||Vu||^2 -/i J diu^di{vl^dku>)dx 


<-f|l|Vn||i2+C||Vu||i3 
<-f|l|Vn||i2+C||Vu||i2||V2u|U2. 
We deduce from integration by parts and (II.Ih that 

h = — VP • {dtu + u ■ Vu)dx 


(3.6) 


/ 


= / PdjuJdiU^dx 


< C\\P\\B,Mo\\c)ju''diU^\\'kii, 


(3.7) 


where one has used the duality of space and BMO one (see [211 Chapter IV]) in the last 
inequality. Since diY{dju) = djdivu = 0 and V-*- • (Vu-^) = 0, Lemma [231 yields 


\h\ = 


PdjuJdiU^ dx 


<C||VP||^2||Vn||i2. 


Then, substituting (|3.6p and (13.81) into (13.5p gives 

f + IIVp«lli= < ^+ Iivp||i2) \\vu\\l2. 

On the other hand, since {p, u, P) satisfies the following Stokes system 

—pAu + VP = —pu, X e 
divM = 0, X e 

^u(x) ^ 0, jxj —>■ 00 , 

applying the standard L^’-estimate to (I3.10p (see [25]) yields that for any r G (l,oo), 

||V^u||Lr + llVPlI^r < C{r)\\pu\\Lr-. 

It thus follows from (13.9p . (13.lip , and (13.11) that 

f |l|Vu||i2 + W^piiWh < C\\^n\\L4^u\\l2 < IwVpuWh + C\\Vu\\\2, 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


which implies that 


J^{p\\^u\\l 2 ) + \\^pu\\l 2 <C\\Vu\\\ 2 . 


(3.12) 
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This combined with (j.S.4l) and Gronwall’s inequality gives (I3.2h . 

Finally, multiplying (I3.12p by t leads to 

I (t/i||Vu||i.) +t\\^u\\h < Ct\\Vu\\l, + C\\Vu\\l„ 

which together with (j3.4p and Gronwall’s inequality yields ()3.3I) and completes the proof of Lemma 

K2[ □ 

Next, motivated by [8lll2lll8j where they deal with the compressible Navier-Stokes equations, 
we have the following estimates on the material derivatives of the velocity which are important 
for the higher order estimates of both the density and the velocity. 

Lemma 3.3 There exists a positive constant C depending only on /r, ||/Oo||LinL°°) IIv^'*^o||l2, and 
II Vuq 112,2 such that for i = 1,2, 


and 


sup f||^ m||^2 + [ f||Vu||^2dt < G, 
te[o,T] Jo 


sup f {\\VW\I 2 + W^PWh) < C. 
tG[o,r] 


(3.13) 

(3.14) 


Proof. First, operating dt + u ■'V to (mii yields that 

dtipiP) + diY{puiP) — p/S.iP = —pdi{diU ■ Vu^) — pdiv{diudiU^) — djdtP — {u ■ V)djP. 
Now, multiplying the above equality by we obtain after integration by parts that 

= — J pdi{diU ■ Vu^)u^dx — J pdiY{diudiU^)vjdx + J 


where 


satisfies 


<G||Vu||i4 + ^||Vu||i2 + J, 


J = — u^dfdjPdx — / ii^u ■ VdjPdx 


J = y PtdiYudx + J u - VPdiviidx + J dju^diPiPdx 
= J {Pt + u- V P)djU^diU^ dx — J Pdju'‘diVp dx 


due to divu = 0. 

Then, since divu = 0, integration by parts implies 
f {Pt + u ■ VP)dju'‘diU^dx 


-f 

dt J 


Pdju'‘diU^dx — I Pdju\diU^dx— / PdjU^diujdx 


/ 


— / Pu^dk[djU^diU^)dx 


d 


= — I PdjU^diU^dx — I Pdju\diU^dx — / Pu^dkdjU^diU^dx 


dt 


/ 


(3.15) 


(3.16) 
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— J PdjU^diuldx — J Pu^dkdiU^djU^dx 

d 


dt 


Pdju'‘diU^dx — / Pdjii^diU^dx + / PdjU^dkU^diU^dx 


/ 


— / PdjU^diiPdx + / PdiU^dkU^djU^dx, 


which together with (|3.16l) . Holder’s and Young’s inequalities yields 


j PdjU^diU^dx + C J |P||Vu||Vu|dx + C' / |P||V 
I Pd,u%u^dx + C {\\P\\l, + ||Vu||i4) + ^llVulli,. 
Next, putting (I3.17p into (|3.15p gives 


J <- 

( 

<- 


ul^dx 


+ | / \Vu\^dx<C\\P\\i, + C\\Vu\\i„ 


where 


satisfies 


^{t) = \Jp\u 


'^dx — / PdnU^djU^dx 


(3.17) 


(3.18) 


- J plupdx — C'||Vu|||^2 <'l'(t) < J p\u\‘^dx + C\\Vu\\\2 

due to (13.8p and (13.lip . Moreover, it follows from Sobolev’s inequality, (13.lip , and (13.11) that 
||H|li4 + ||Vu||i4 < C (||VP||t4/3 + 

< c\\pu\\%,, < cMU^uWl^ < cw^iiWh. 


(3.19) 


(3.20) 


Multiplying (|3.18p by P{i = 1,2) and using (|3.19P " (|3.2Up . we obtain (I3.13P from Gronwall’s 
inequality, (13.21) . and (13.3p . 

Finally, (I3.14p is a direct consequence of (I3.13P and (13.IIP . The proof of Lemma l3.3l is finished. 

□ 


3.2 Higher Order Estimates 

The following spatial weighted estimate on the density plays an important role in deriving the 
bounds on the higher order derivatives of the solutions {p,u,P). 

Lemma 3.4 There exists a positive constant C depending only on p, ||/Oo||Lini/°°; IIpo^“||l 1; 
||y^Mo||L 2 , ||Vuo||x, 2 , Nq, andT such that 

sup ||px“||ii < C. (3.21) 

i6[0,T] 

Proof. First, for > 1, let y ?at € C'^(M^) satisfy 

0<7PJV<1, (/p^(x) = J \VipN\<CN-\ (3.22) 

10, |x| > N, 

It follows from (ll.ip i that 


fL 

dt 


ppi^dx = / pu-Vifi^dx 
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> -CN 


-1 


X 1/2 / . X 1/2 

pdx\ ( / p\u\'^dx\ '>—CN~^, 


(3.23) 


where in the last inequality one has used (I3I]) and (133]). Integrating (I3.23p and choosing N = 
Ni = 2 Nq + ACT, we obtain after using (|1.4p that 


inf / pdx > inf / ptpN.dx 
0<t<Tj 


> J poifiNidx — CN^^T 

> [ podx - , 

2No + ACT 

> 1/4. 


(3.24) 


Hence, it follows from (I3.24p . (13.ip . (12.2p . (13.4h . and (|3.2I) that for any p G (0,1] and any s > 2, 

Wux-^L^/r, < C (||p^/\||l2 + ||Vnb2) < C. (3.25) 

Multiplying mi by x°' and integrating the resulting equality by parts over yield that 


_d 

dt 


< CWpX^" 8+0 II s+a\\ux 8 +“||r 8 +a 
— nr "iTTs" 


<C J px'^dx + C, 

which along with Gronwall’s inequality gives (I3.2ip and finishes the proof of Lemma 13.41 
Lemma 3.5 There exists a positive constant C depending on T such that 

sup ||p|binivi..j + / fl|V^M||i2 + ||V\|U +t\\V‘^u\\l2^^^'] dt 
te[o,T] Jo \ J 

+ £ (llVPlli^ + ||VP||^ + dt < C(T). 

Proof. First, it follows from the mass equation mi that Vp satisfies for any r >2, 

d 


□ 


dt 


\\Vp\\Lr<Cir)\\Vu\\L^\\Vp\\Lr. 


(3.26) 


(3.27) 


Next, one gets from Gagliardo-Nirenberg inequality, (|3.2p . and (13.lip that for g > 2 as in Theorem 

m 

0 — 2 a a 

(3.28) 


0-2 Q 

I 2(0-1) ||V72„,|| 2(g-l) 


liVnIUco < C||Vn||2r^M|V"u||2r'' < c||ph||^. 

It follows from ()3.24p . (13.ip . (12.21) . and (|3.21l) that for any p G (0,1] and any s > 2, 

Srja SrjcL 

Wp'^vW s <C\\p'^X^\\ 4s I lux 4V|| 4s 
Wj^n — lie 11^3^ II iij^— 

< C||H|)5^I|4 S“I|S (||p‘'"o|Il» + llVolli. 


< c (||P+t.||i,. + ||Vi,||i«) , 


(3.29) 
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which together with the Gagliardo-Nirenberg inequality shows that 


2(9-1) 

ij2_2 


\\pu\\l<i < C\\pu\\^2 \\p 


9 ( 9 - 2 ) 
I 9^-2 


2(9-1) 


< C\\pu\\^2 (|lv^^i||L2 -h \\Vu\\l2) 


9(9-2) 

2-2 


2(9-1) 

92-2 


<C + llVull/, 


9(9-2) 

92-2 


(3.30) 


This combined with (j3.2l) and (I3.13li leads to 


rT 


9+1 


IIp^IIl' +4P^L‘i]dt 


q'=’-Q^- 2 q-l 


< C* y ^11 VP^IIl2 + t||Vh||^2 + t -pljdt 

<C, 

which along with (|3.28h in particular implies 

[ llVullioodt < C. 

Jo 

Thus, applying Gronwall’s inequality to (|3.27l) gives 

sup ||Vp||i2nL9 < C. 

4g[0,T] 

Furthermore, it is easy to deduce from (13.lip . (I3.3ip . (13.ip . (13.2h . and ()3.3p that 

rT / 9+1 \ 

|V^tt|||2 + + t||V^u|||2ni9 ) dt 
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+ 


g+1 

|VP||i 2 + ||VP||^ 9 ^ + t||VP||i2nL9 ]dt<c, 


(3.31) 


(3.32) 


(3.33) 


(3.34) 


which together with (|3.1I) and (|3.33l) yields (|3.26l) and completes the proof of Lemma 13.51 □ 

Lemma 3.6 There exists a positive eonstant C depending on T such that 

sup Wpx'^WL^nH^nw^’i < C{T). (3.35) 

te[o,r] 

Proof. One derives from pi.ip i that satisfies 

dt{px°‘) + u ■ V(px“) — apx°‘u ■ Vlogx = 0. (3.36) 

Taking the Xj-derivative on the both side of (13.361) gives 


0 =dtdi{px°‘) + u ■ Vdi{px°‘) + diU ■ V{px°‘) — adi{px°‘)u ■ Vlogx 
— apx°‘diU ■ V log X — apx°‘u ■ diV log x. 


(3.37) 


For any r G [2, g], multiplying (I3.37p by |V(px“)|'’ ‘^di{px°') and integrating the resulting equality 
over we obtain after integration by parts that 

^||V(px“)||l. <C (1 + llVulUoo + ||U • Vlogxllioo) ||V(px“)|U. 


10 



























+ CII/ox^IIloc (|||Vti||Vlogxl||Lr- + Ill'll Ilog xlllir-) 

<C{l + \\Vu\\wi.)\\V{px^)\\Lr 
+ C 11/9X^11^00 


_ 2 ||__ 3 .. 

Lr- + \\UX 5 \\j^4r\\X 2 || 4 ^ 

L ^ 


<C (1 + ||V?x||vi/i.9) (1 + ||V(px“)||l’- + ||V(px“)||l9) , 


(3.38) 


where in the second and the last inequalities, one has used ()3.25p and (I3.2ip . respectively. Choos¬ 
ing r = g in p3.38D . together with (|3.26D . indicates that 

sup ||V(px“)||l9 < C. (3.39) 

iG[o,r] 

Setting r = 2 in (|3.38D . we deduce from (13.2611 and (13.3911 that 

sup ||V(/9x“)||i2 < C. (3.40) 

iG[0,r] 

This combined with (I3.39p and (I3.2ip gives (|3.35ll and finishes the proof of Lemma 13.61 □ 

Lemma 3.7 There exists a positive constant C depending on T such that 

sup t\\y/put\\l 2 + [ t\\Vut\\l 2 dt <C{T). (3.41) 

te[o,T] Jo 

Proof. Differentiating 0)2 with respect to t leads to 

putt + pu ■ Vut - pAut + VPt = -ptut - {pu)t ■ Vu. (3.42) 

Multiplying (I3.42p by ut and integrating the resulting equation by parts over we obtain after 
using Oi and (jl.ip -? that 


2 dt 


\\^/put\\‘l 2 + p\\'^ut\\‘l 2 = - j pt\ut\'^dx - j{pu)t ■ Vtt • utdx 

< Cy p|u||ttf|(|Vrtt|-|-|Vtt|^-|-|tt||V^ti|)(ix 

+ C y p\u\‘^\Vu\\Vut\dx + C y p\ut\‘^\Vu\da 

AK 1 + K 2 + K 3 . 


(3.43) 


We estimate each term on the right-hand side of (I3.43P as follows. 

First, the combination of (I3.25p . p3.29p . (13.4p . and (13.2p gives that for any p £ (0,1] and any 
s > 2, 


Wp'^nWis/r, + \\ux "’Wls/v < C, 


(3.44) 


which together with p3.29p . p3.2p . and Holder’s inequality yields that 

K, < C\\^uhe\\^Ut\\]!2"\\VpM]!e {W^Uth^ + WVuW^) 

+c\\p^/\\\it2\\vput\\];,^\\4/^^^^ 

< C\\^ut\\]!2" (WVpMl^ + l|V«i||L2)i/2 (||V«,||^, + \\vMl^) 

< ^ll^'*^t|lL2 -I- (7 (1 -I- \\y/put\\'^^2 + ||V^tt||^2) . 

Next, Holder’s inequality combined with p3.29p . (j3.2p . and (I3.44p leads to 


(3.45) 


K 2 +K 3 < C\\y/pu\\ls\\Vu\\L4\\Vut\\L2 + W'^uWl^WVpMl^WVpMl^ 


, 1/2 
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(3.46) 


^ + C* (l + \\y/put\\\2 + ||V^u|||2) • 


Substituting (|3.45l) - (l3.46p into (j3.43p gives 
d 


dt 


W^putWi^ + p\\yut\\i 2 < C\\^put\\i 2 + C\\V^u\\i 2 + C. 


(3.47) 


Finally, noticing that it follows from (13.21) . (I3.44p . (13.261) . and Garliardo-Nirenberg inequality 


that 


[ \\^/put\\hdt< [ {\\^u\\l2 + \\^/p\u\\Vu\\\l2)dt 
Jo Jo 

<C + C [ \\,/^u\\l,\\Vu\\l,dt 

Jo 

'[ {\\Vu\\l 2 + \\vMh)dt 
Jo 


<C + C 


<C, 


(3.48) 


multiplying (13.471) by t, and using Gronwall’s inequality, (I3.26p . and (I3.48p . we obtain (13.411) and 
complete the proof of Lemma 13.71 □ 


4 Proof of Theorem 11.11 

With the a priori estimates in Section [3] at hand, we are now in a position to prove Theorem ll.il 
By Lemma 2.1, there exists a T* > 0 such that the problem (ll.ip - (jl.2p has a unique local 
strong solution {p,u,P) on x (0,T*]. We plan to extend the local solution to all time. 

Set 

T* = sup{T I {p,u,P) is a strong solution on x (0, T]}. (4.1) 

First, for any 0<r<T*<T<r* with T finite, one deduces from (13.21) . (I3.14p . and (13.411) that 
for any q > 2, 

Vu e C([r,r];L2nL''), (4.2) 

where one has used the standard embedding 

L°^{T,T;H^)nH\T,T;H-^) ^ C{t,T-,L^) for any q£[2,oo). 

Moreover, it follows from (I3.26p . (I3.35p . and [T9l Lemma 2.3] that 

p£C{[0,T]-,L^nH^nW^’'^). (4.3) 

Finally, if T* < oo, it follows from (14.21) . (|4.3I) . (13.2p . p3.4l) . and (I3.35p that 

{p,u){x,T*) = lim {p,u){x,t) 
t^T* 

satishes the initial conditions dLSI) at t = T*. Thus, taking {p,u){x,T*) as the initial data. 
Lemma 2.1 implies that one can extend the strong solutions beyond T*. This contradicts the 
assumption of T* in (14.ip . The proof of Theorem 1 1.1 1 is completed. □ 
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